Perturbative expansion for the half-integer rectihnear 
dischnation hne in the Landau-de Gennes theory 



H. Arodz* and R. Pelka^ 
February 2, 2008 



Abstract 

The structure of the half-integer rectihnear dischnation hne within the framework of 
the Landau-de Gennes effective theory of nematic hquid crystals is investigated. The 
consistent perturbative expansion is constructed for the case of L2 7^ 0. It turns out 
that such expansion can be performed around only a discrete subset of an infinite set of 
the degenerate zeroth order solutions. These solutions correspond to the positive and 
negative wedge dischnation lines and to four configurations of the twist dischnation 
line. The first order corrections to both the order parameter field as well as the free 
energy of the dischnation lines have been found. The results for the free energy are 
compared with the ones obtained in the Frank-Oseen-Zocher director description. 
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1 Introduction 



Disclination lines in nematic liquid crystals rank among the most popular topological defects 
encountered in condensed matter physics P^ISEllll- This notwithstanding, there are still 
a number of interesting theoretical problems to be solved, which particularly holds for the 
Landau-de Gennes theory, for its review and discussion see, e.g., [Sj. Such a state of affairs is 
partly due to the relative complexity of the order parameter structure in that theory. Not less 
importantly, the Landau-de Gennes model appeared a decade later than the Frank-Oseen- 
Zocher director formalism In contrast to that formalism, which leaves the question of the 
core structure of the disclination line open, the relevant solution in the Landau-de Gennes 
model is smooth everywhere [3 El El Although both the formalisms account for the 
intrinsic elastic anisotropy of a nematic medium, only in the director formalism the exact 
solution for generic values of the Frank elastic constants had been found fUlEI- At the same 
time, to the best of our knowledge, even approximate analytical results for the case of L2 7^ 
in the Landau-de Gennes model, which correspond to the two-constant approximation in 
the director formalism, are still missing. The only results are those reported by Schopohl 
and Sluckin in jH] , where numerical calculations have been carried out for half-integer wedge 
disclinations. 

In this paper we consider a generic half-integer disclination line within the framework 
of the Landau-de Gennes theory. Our main result is the formulation of a consistent pertur- 
bative expansion for the case of L2 7^ 0. It must be mentioned that we draw extensively 
on the results obtained in ^0]- Therein the isotropic L2 = case has been thoroughly dis- 
cussed, which represents the zeroth order of the perturbative scheme constructed here. Via 
some symmetry considerations we find that only around a restricted set of the zeroth order 
solutions a perturbative expansion can be consistently performed. The presented scheme 
allows one to generate the corrections in all consecutive orders of the expansion. The forms 
of the first order contributions to the order parameter have been found explicitly. For the 
positive and negative wedge disclination lines in the simplest high symmetry case of /5 = | 
accounted for in [lOj also numerical calculations have been performed. What is important, 
the relevant corrections proved to be at least one order of magnitude smaller than the zeroth 
order contribution, which validates the apphcability of the preturbative approach. Next, 
the polynomial approximation suggested in [THj has been taken full advantage of in order to 
estimate the first order corrections to the free energy. In accord with the results obtained 
within the two-constant approximation in the director formalism [T2], for L2 > (L2 < 0) 
the twist disclination line (the positive wedge disclination line) has been found to have the 
lowest free energy. The first order correction to the free energy of the twist disclination 
line was found to be logarithmically divergent, whereas that of both the wedge disclination 
lines turned out to be finite. The corresponding first order corrections to the value of the 
'critical' temperature below which the smooth disclination line is energetically favourable 
have been estimated. The corresponding 'critical' temperatures have been found to increase 
(decrease) for L2 > (L2 < 0) for all types of the dischnation line. 

The paper is organized as follows: in Sec|2] the preliminary considerations have been 
presented. Sec 01 includes the very construction of the perturbative expansion, whereas SecH) 
focuses on the first order corrections to the free energy and to the value of the 'critical' 
temperature, in SecIHjsome final remarks have been made, and in Appendix the derivations 
of the first order corrections to the order parameter are to be found. 
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2 Preliminaries 



2.1 The Landau-de Gennes theory 

The order parameter for the nematic hquid crystal in the Landau-de Gennes theory is a 
symmetric traceless real tensor Qij, where i,j = 1,2,3 and the corresponding free energy 
density is given by the following formulae 

:F = ^LidiQjkd.Qjk + ^L2diQ,kdjQjk + V{Q), (1) 

where 

V{Q) = -|Tr (g^) - ^Tr (Q^) + | (Tr {Q^))' . (2) 

For the standard stability reasons constants a and c are assumed to be positive. Fur- 
thermore, for most of the explored nematogenic compounds, see e.g. I3j, the ground state 
proved to be that of a uniaxial ordering, that is why we set also the constant b to be positive. 
Simple stability analysis for basic inhomogeneous conformations of the order parameter jT3] 
shows that the elastic constant Li should be positive, whereas for L2 the following restriction 
has to be satisfied 

L2 > -h,. (3) 

Hence, negative values of that constant are not excluded. 

The order parameter in the ground state is a constant tensor Qg, whose form we obtain 
from the necessary condition for local extrema of the free energy. It is given by the following 
formula 



Qg = OQoO' 



(4) 



where 



Qo = Vo 




(5) 



and 



Vo 



b + Vb'^ + 24ac 



(6) 



whereas the transformations O without loss of generality can be restricted to the set of 
the proper three-dimensional rotations 5*0(3). Thus the vacuum manifold is continuously 
degenerate. 

In it has been found that for the rectilinear disclination line in the Landau-de Gennes 
theory with L2 = the order parameter takes on the following axially symmetric form 



= -VoO{<P) 



/ ^ + Rip) 

^ - Rip) 

V 








o^i^), 



(7) 



■ISip) 



3 



where tiq is given by the formula functions S{p),R{p) are a pair of smooth structure 
functions satisfying the following boundary conditions 



S{0) = const, S{oo) = 1, (8) 

R{0) = 0, R{oo) = 1, (9) 

and the matrix 0{(j)) has the form 



/ cos f - sin f 



sin I cos f I . (10) 
V 



The core of the disclination line given by the above solution is smooth with a planar uniaxial 
phase on the line, the maximal eigenvalue of Q being negative, reached via the ring of 
biaxiality as one approaches the line from the far-off area filled with the uniaxial phase. On 
taking some simple assumptions as to the creation of the singular disclination core it 
was pointed out that for sufficiently low temperatures parametrized by 

the smooth disclination line represents an energetically preferred conformation, wheras for 
/3 large enough a singular core with the isotropic phase inside will be more favourable. The 
'critical' value of (3 dividing those two stability domains was estimated to be equal to 0.12. 

2.2 Degeneracy of the ground state of a rectilinear disclination 
line 

It is clear that Ansatz (jZj) does not exhaust the set of the possible solutions of the pertinent 
Euler-Lagrange equation for the rectilinear m = | dislination line when L2 = 0. Other 
topologically equivalent solutions are given by the formula 

Q{p,^) = OMp,<P)Ol, (12) 

where d G SO (3). 

The topological equivalence does not yet mean the physical one. While all the solutions 
fll2|) are characterized by the same free energy, they describe distinguishable spatial config- 
urations of the nematogenic molecules; e.g. Qd corresponds to a configuration where far off 
the disclination line all the molecules lie on avarage in the xy-plane (a wedge disclination 
line), whereas inserting for Oi in formula (|T2|l the matrix 




(13) 

one arrives at the configuration referred to in the literature [2j as a twist dislination line 
where all the far-off molecules tend to lie in the plane parallel to the direction of the 
line. On the other hand, among solutions ()12|) one may find such which, failing to display 
mathematical identity, result in the same spatial configurations of molecules, the difference 
originating from the mere change of the frame of reference. In general, two solutions Q and 
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Q give physically equivalent spatial configurations, if there exists such a global coordinate 
transformation Oc € Hd that 

= Q- (14) 

Group Hd contains all proper orthogonal transformations which leave the position of the 
disclination line unchanged. We purposefully omit the group of one-dimensional transla- 
tions along the line as irrelevant here. Thus Hd = S0{2). Eq. ()14|1 means that physically 
nonequivalent configurations correspond to the different right cosets S0{3)/ S0{2). Know- 
ing that any rotation can be parametrized with the set of three Euler angles, one finds easily 
that the looked for elements of the coset are given by the formula 



Oi(e,^) = o,(e)a(^), (15) 

where Oy, Oz refer to the rotations about the y- and z-axis, respectively, and 9 e [0, tt], e 
[0,27r). One should also exclude the case where both the rotation matrices in the above 
formula commute, i.e. when ^ = or tt, and different values of angle \l/ lead to the physically 
indistinguishable configurations. Apart from that, there exists one more symmetry in the 
system. Namely, one can easily check that 

Qd = OoQdO^, (16) 

where 



Oo=\ -10 . (17) 




This means that transformations Oi{Q, \1/) i Oi(6, '^)Oo result in mathematically equivalent 
solutions. Moreover, because 



0,{e,^)0o = 0i{e,n + ^), (18) 

we get nonequivalent solutions for \E' in the reduced segment of [0,7r). Summing it all up, 
the physically nonequivalent configurations are given by formula p2|) . where Oi is given 
by formula ()15|1 with pairs (0, \I') lying in the interior of a square (0, vr) x (0, vr) or on its 
boundary segment [0, vr] x {0}. 



3 Construction of the perturbative expansion 
3.1 First steps 

The main goal of this section consists in reporting the construction of the perturbative 
scheme in the case of L2 7^ 0. The unperturbed system corresponds to the so called one- 
constant approximation in the director formalism, where Ki = K2 = = 9?7g(2Li + L2), 
whereas the perturbation takes on the following form 

(5^ = ((divg)2 - ^(gradQ)^) , (19) 
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where (divQ)^ = diQijdkQkj, and (gradQ)^ = dkQjidkQji- As the expansion parameter we 
identify the dimensionless ratio e = 2L2/ (2Li + L2). One can easily check that the following 
identity holds: 

2(K-K2) , , 

e = ^ ^ '\ (20) 

where K = {Ki+K^) /2. Using the experimental data TT for two popular representatives of 
nematic liquid crystals (for MBBA at transition temperature ~ 320 K , Ki = 7 .1^ 10^^^ A^, 
K2 = 4.0 X 10-^2 N, = 9.2 X 10-^2 jy. for DIBAB at transition temperature ^ 305 K, 
Ki = 4.8 X 10-12 = 3.0 X lO'^^ jy^ ^ 4 7 ^ lO'^^ arrive at the following 

estimations 

euBBA ~ 1-01, ^DIBAB ~ 0.74. 

The fact that these values of e are large does not mean that the perturbative calculations 
are out of place here, as one should remember that what is relevant is not the value of the 
expansion parameter itself but the order of magnitude of the full consecutive corrections. It 
will be demonstrated later that the first order corrections are at least one order of magnitude 
smaller than the zeroth order contributions, so the perturbative approach finds quite a good 
justification. 

For the sake of operational convenience let us introduce a rescaled order parameter q 
defined by the following formula 

Q = Ivoq, (21) 
as well as dimensionless coordinates defined as follows 



= ^os\ z = 1, 2, 3, gdzie = • (22) 

In the variables introduced above the Euler-Lagrange equation takes on the following form 



dkdkQij - le \2{dkdiqki)Sij - 3 (didkqkj + djdkqki) + Sdkdkqij] + ^qij + 3(3{q 



-i(l + 3/3)Tr(g2)g^^. _ fiTr{ct)b,, = 0, ' 



(23) 



where the derivatives with respect to the dimensionless coordinates are marked out by a 
tilde. We seek its solutions assuming the perturbative expansion 

q = qo + eqi+ e^2 + ■■■■ (24) 

Inserting expansion into equation (j^ . and equating to zero the factors appearing 
at the subsequent powers of e we find equations governing the dynamics of the subsequent 
corrections to the order parameter. 

In the zeroth order we obtain the equation discussed in detail in [TU] 

dkdkqo + 1% + mo - ^(1 + 3/?)Tr(g2)go - /3Tr(g2)/ = 0. (25) 

According to what has been said in the previous section the physically nonequivalent 
solutions of equation are parametrized by a pair of angles (0,\1') with values in 
(0, vr) X (0, tt) U [0, tt] x {0}. They have the form 
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= 0^{Q,^)Us,4>)0UQ,^), (26) 
where Oi{Q, \&) is given by the formula (fTH|l and qd{s, 0) has the following form 



/ is (s) + R{s) cos ^ R{s) sin (p 
i?(s)sin0 I S (s) — R{s) cos (j) 

V -ls{s) 



(27) 



The equations in higher orders of the perturbative expansion are of the general form 



{2i 



where k = 1,2, . . ., C denotes a linear operator in the space Ai of coordinate-dependent 
symmetric traceless and real matrices 3x3. Its action on a generic element m = m(s*) in 
that space is defined by the formula 

~ ~ 2 1 1 

Cm = dkdkrh+ -m+?,[3{%m+m^ 

(29) 

The consecutive terms on the r.h.s. of Eq. (j^H|) can be obtained from the formal recipe 

1 d'^h 



where matrix h is defined by the formula 



k\ de^ 



(30) 



e=0 



hij{s'] e) = le[2{dkdkqki)Sij - 3{didkqkj + djdkqki) + ^dkhqij] - 3P{f - m - q%)ij+ 
+i(l + 3/?) [Tr(g2)g - 2Tr{m)qo - Tr{ql)q]i, + /?Tr[g(g - 2q,)]5,,, 



with q in the form of expansion For example 



(31) 



= - [2{dkdiqoki)Sij - 3{didkqokj + djdkqoki) + ^dkdkqoij] 

= ^[2{dkdiqiki)Sij - 3{didkqikj + djdkqiki) + ^dkhquj] - 3(3{ql)ij 

+ -{1 + 313) \Tr{ql)% + 2Tr(gogi)gil . . + l3Tr{ql)5i,. 



(32) 



(33) 



C and A/fc depend on the zeroth order solution (jHIIEU)- In spite of that, we do not fit 
those entities with an additional subscript referring to a given go^'*^ but stipulate that all 
the formulae to be found hereafter refer to a generic choice of a pair (B, \E'), unless stated 
otherwise. 



3.2 Symmetry considerations 

It is clear that the unperturbed system is invariant under the action of the full 50(3) group, 
i. e. the free energy is the same for any solution of the form 



go = Oq,o^ 



(34) 
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where O G 5*0 (3). A simple consequence of this symmetry is the existence of the corre- 
sponding rotational zero modes. One can easily check that they have the form 



Oa = [ia,qo], (35) 

where ta denotes the generator of the rotation about axis a {a = 1, 2, 3). The modes (jH^j) 
are eigenfunctions of operator C belonging to the eigenvalue equal to zero. 

Yet, the 50(3) symmetry of internal rotations is not the only one in the system. The 
other one is the group of translations along the three independent directions. This symmetry 
originates from the homogeneity of the system, as we assume here that all the material 
constants a,b,c, Li, L2 do not depend on the coordinates x\ Similarly, there exist three 
corresponding translational zero modes 

= (36) 

where i = 1,2,3. Furthermore, because go does not depend on variable s^, the mode fs 
corresponding to the trivial translations along the disclination line vanishes. On differenti- 
ating both sides of equation (j2Sl) with respect to {i = 1, 2) we see that these modes are 
the eigenfunctions of C belonging to the eigenvalue 0. 

The zero modes discussed above possess a clear physical interpretation, i. e. they appear 
as the consequence of the symmetry of the unperturbed model. 

Note that 

[A4]=0 i 4^ = 0. (37) 

From that follows that there exists such a solution of equation ()28|1 with k = 1 that is inde- 
pendent of coordinate s^. However, the general solution of that equation does not exclude 
the s^-dependence. That is, a particular solution dispalying the translational symmetry 
along the direction of the line may in general be appended by a solution breaking that 
symmetry. The latter one would have to be an eigenfunction of C to the eigenvalue and 
would correspond to the excitations of the disclination line disturbing its rectilinearity. As 
the discussion of such perturbations will clearly exceed the intended scope of the present 
paper, we set that part of the full solution to zero in all consecutive orders of the perturba- 
tive expansion. Now the solution depends exclusivly on coordinates s^, s^. So, everywhere 
hereupon the Roman indices take on only two values {1,2}. 

Now that we have explored the symmetries of the unperturbed model we may pass on to 
discussing their consequences. Let fi denote any of the zero modes mentioned above. Next, 
multiplying both sides of Eq. ()28|) by fi, taking the trace, and integrating the resulting 
equation over the circular area C{0,L) centered at the origin of the coordinate system, 
whose finite radius we denote by L, we arrive at the following equation 



d^sTiifiCqk) = f (fsTiifiAfk)- (38) 

/C(0,L) Jc{0,L) 

The l.h.s. of the above equation can be rewritten into the form 



l.h.s. = / d^s Trii^idAqk) - {didifi)qk + qktfi]. (39) 

JC(0,L) 

The contribution from the last term in brackets vanishes, as (x represents a zero mode, 
whereas the contribution from the remaining terms, on applying the Stokes theorem and 
some simple algebra, may be transformed into the following form 
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where s 



-A\2 



+ S 



2^2 



Here we arrive at the crucial point of this paper. Firstly, note 
that the rotational zero modes da reach their 0-dependent asymptotic form as quick as 
whereas the translational ones vanish as s~^, see formulae (j2Z|), Secondly, it is only 

natural to assume that the consecutive corrections are finite at infinity. That given, we 
conclude that the expression pn]) vanishes in the limit L oo. As a consequence we obtain 
the following consistency conditions 



lim / d^sTiiflAfk) = 0, 

L^oo Jc{0,L) 



(41) 



where k = 1,2, . . .. They are to play an essential role in the construction. 

Let us analyse the consistency conditions (pT|l for k = 1. To this end it is convenient to 
write J\fi in the form 
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2Tr(Do)/ - 3(1^0 + D^) + 3Ago 



(42) 



where Do denotes the following matrix 



{Do)^, = dSi {qo)ij , (43) 

and A is the Laplace operator in the plane {s^, s^). Remembering that the zero modes are 
symmetric and traceless matrices, the consistency condition can be rewritten as follows 



lim / d'^sTi 

L^oo Jc{0,L) 



fi (2Do - Ago 



0. 



(44) 



To begin with, let us consider the translational zero modes (i = 1,2). Integrating by 
parts one can transform the l.h.s. of the condition into the following form 



lim <^ / [26kn{ri)ji{rn)kj 

L^OD [JdC{0,L) 



d s di 



C{0,L) 



[ri)jk[ri)kj 



dd' 



di{qo)ijdk{qo)kj - ^dk{qo)jidk{qo)ij 



(45) 



As the expression within the brackets in the second integral vanishes in the limit s ^ oo, 
this integral will give no contribution. Neither will the surface integral, which follows 
immediately from the asymptotic behaviour of fj's. Thus, the consistency conditions with 
the translational modes are trivial; they do not impose any restrictions on the zeroth order 
solution. There is nothing strange about it because the unperturbed model as well as the 
perturbation itself do not break the translational symmetry, so this symmetry holds for the 
full system. As a consequence, we obtain an infinite set of degenerate solutions which differ 
from each other in the position of the disclination core. 

The consistency conditions with the rotational zero modes 9a (a = 1, 2, 3), see Eqs. (jHSl 
HT|) . can be rewritten as follows 



lim / d^sTi{ia(2[qo,Do]-[qo,Aqo])} =0, 



(46) 
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where the brackets denote the commutator. One can easily see that [go, ^Qo] = 0. For the 
remaining term, on inserting the exphcit form of q (Eq. (j26|) ) and using polar coordinates 
in the plane {s^, s^), we obtain 



lim |--sin^esin(4*) [^dssR\R" + 

L^oo [2 Jo \ 



R!_ 

s 



R 



lim < sin G cos G / dss 

L->oo [3 JO 

= 0, 



sis" + -S' 

s 



R' R' 



3RiR" + 



(47) 

(48) 
(49) 



for a = 1,2,3, respectively, where ' denotes the differentiation with respect to s. The 
integrals in (jTfjl and P8|l are nonzero and finite for every finite radius L, whereas in the 
limit L ^ oo they diverge. Therefore, the only way to satisfy the consistency conditions is 
to select these pairs (G, \E') that solve the system of equations 



sin^ G sin(4\E') 
sin G cos G = 



0, 



(50) 



which has the following solutions 



(51) 
(52) 

(53) 
(54) 
(55) 
(56) 

Thus, via the consistency conditions we have arrived at the conclusion that not all configu- 
rations are in agreement with the perturbation. The consistency conditions serve as a kind 
of a selection rule for the zeroth order configurations 'matching' the specific perturbation 
()19|). The situation here is similar to what one encounters while considering perturbations 
of the quantum mechanical systems with a degenerate level, where the perturbation is well 
known to single out a basis in the subspace of the Hilbert space corresponding to that level. 
So, each of the pairs (G, given above represents the zeroth order configuration singled 



G = 


0, 




= 0, 


G = 






= 0, 


G = 


TT 

2' 




= 0, 


G = 


71 

2' 




TT 

" 4' 


G = 


TT 

2' 




TT 

" 2' 


G = 


TT 

2' 




3tt 
~ T 



out by the perturbation 6 J-". The first two configurations ()51l52j) correspond to the wedge 
disclination lines, positive and negative, respectively, whereas the last four cases are twist 
configurations. 

Let us comment on the result for a = 3, see Eq. (jinjl . That identity suggests that 
the subgroup SO (2) of the rotations about the direction of the dischnation line may well 
represent the symmetry of the full system. Yet, one can check that, in general, those 
rotations do affect the value of the perturbation. So, all one can gather from Eq. ()49p is 
that they do not change the free energy up to the first order of the perturbative expansion. 

In conclusion, it is only some of the zeroth order solutions ql^'^^ that are 'matched' 
with the considered perturbation. In other words, the perturbative expansion of the or- 
der parameter will be consistent only when carried out around the selected configurations 
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given by ()51II56|) . In the next subsection we present a systematic procedure for finding the 
consecutive corrections to the order parameter. 



3.3 Corrections to the order parameter 

The corrections to the order parameter field in the consecutive orders of the perturbative 
expansion represent the solutions of the linear differential equations ()28|) . where the oper- 
ator £, given by is the same for a fixed go in all orders of the expansion, and the 
inhomogeneous terms A4 {k = 1,2, . . .) are defined in (fHUjl . (jHH). 

Let us begin with the construction of the solution for 6 = \E' = (a positive wedge 
disclination line), and = vr, \E' = (a negative wedge disclination line). Next we will 
pass on to discussing the remaining cases ()53II56|) . It is convenient to introduce the local 
0-dependent basis in the space Ai defined by the following series of formulae 



£;:(0) = O(0;a)E°O^(0;a), (57) 



where a = 1, . . . , 5, and 



cos I —0" sin I 



O(0;a) = I ffsinf cos f | , (58) 
1 



while 



s;o _ j_ 












r 1 









1 


i 


EO - j_ 





-1 


") 


\ 







^ 





/ 



-C/4 



1 

V2 



1 ° 















:) 


E^ - ^\ 


V 1 









(59) 



(60) 



The index a = ±1 corresponds to the case of the positive and negative wedge disclination 
line, respectively. The set {E^{(f)), a = 1, . . . , 5} is an orthonormal basis with respect to the 
standard scalar product in matrix spaces 



(Qi|Q2) = Tr(QiQ2). 



(61) 



We find it convenient to adopt here and in what follows the quantum mechanical bra-ket 
notation. 

Next step consists in projecting Eq. (j^H|l onto the directions of the basis vectors, which 
results in the following set of equation 

= (KlAT^), (62) 

where a = 1, . . . ,5, k = 1,2, . . ., and there is no summation with respect to a. Inserting 
between and qk the appropriate identity operator \E^){E^\, where the summation is 
carried out exclusivly over b = 1, . . . , 5, one obtains a system of five equations 



5, 



(63) 
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where 



'-'ah 



. = (km = M<ikk), 



(64) 



(65) 



It is clear, see fl29j) . that the matrix elements of are not simple scalars but they contain 
differentiation with respect to polar coordinates (s,0). Some careful algebra allows one to 
obtain 



^Ib = Sat As + . 



1 



1 2 9 



(66) 



where a,b = 1, . . . , 5, 5ab is the Kronecker delta, = ("S^)) and matrices Ji, J2 have 
the following form 



(67) 





/ 





\ 




/ 








\ 




1 










0-1 















1 





, J2 — 


1 















\ 













1 

2 












\) 







1 
2 





/ 


Non-vanishing matrix elements 


of J3 


are 


given by the formulae 









311 



^-2/35- 1(1 + 3/5) (S2 + i?2 
J312 = J321 = 75i?[6/3 - (1 + 3/3)5] 
1/322 
1/333 

^344 = I - /?(5 - 3i?) - |(1 + 3/3) (52 ^ 3^2 

^355 = I - /?(5 + 3i?) - 1(1 + 3/?)(S2 + 3i?2 



2 + 2/35-i(l + 3/3)(52 + 9i?2), 
+ 2/3^ -I (1 + 3/3) (52 + 3/22), 



(68) 



where parameter /3 is defined in (jTTjl , and functions 5", R define the zeroth order structure 
of the order parameter, see ((Zj). 

Note that the only ^-dependence of the matrix elements L"^^ in the basis E"^ lies in 
differential operators That leads to a great simplification in the next step, where 

the Fourier analysis of components Qkai^^^) ^^^^ t)e undertaken. Moreover, the formulae 
()66ll68p show that the full space 3 splits into the simple sum of two subspaces 
A4 = Ai^^^ (B A^*^2) invariant with respect to the action of therein. The subspace Ai^^^ 
comprises components q^^ with a = 1,2,3, whereas A^'^^) _ i]^q remaining ones. We will 
address the subspaces separately. 

For each k = 1,2, 
series with respect to 



we expand components {q^a^^^ = 1)2,3} of A4^^^ in the Fourier 



Qkais, 0) = <(s) + E [<^(s) cosm + <i(s) sin(/0) 



(69) 



1=1 



The Fourier analysis results in a fine structure of an infinite-dimensional space which 
splits into an infinite simple sum of subspaces J-'^^^^^ invariant with respect to the action of 
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C'^^^\ where superscript (1) refers to that part of which acts in One may write 

then 



1=0 



(70) 



The subspace JFq has the following structure 



J^^^^ = \SEl 



(71) 



where S denotes the space of smooth functions of the radial variable s, whereas the structure 
of the subspaces with I ^ can be described as 



■<j{i) 



SE^cos{l(p) © SE^cos{l(j)) © SElcos{l(j)) 
SElsiTL{l(t)) \ © \SE^ sin(Z0)l © \SEl sin(/0) 



(72) 



So, the three-dimensional subspace Tq'^^^ comprises components {vl^{s), a = 1, 2, 3}, whereas 
the components belonging to the six-dimensional subspaces J^^^^ (/ ^ 0) are {v'ka{s),u'{}^{s), 
a — 1,2, 3}. The action of JC,'^^^^ in J^^^^ is represented by which for Z 7^ are 6x6 

matrices with the following block structure 

/ 



^a(l) _ 
H — 



(73) 



where / in the formula above denotes a 3 x 3 identity matrix, whereas , and 
represent submatrices of J\, J2 and J^, respectively, comprising the first three rows and 
columns. For Z = we obtain the 3x3 matrix given by the following formula 



1 



(1) 



(74) 



Note that this operator is independent of cr. 

The last step is a partial diagonalization of C'^^^'^ with Z 7^ 0. It is carried out so as to 



get rid of the off-diagonal blocks containing J2'\ Because C'^^''' are symmetric {[J2 j 

—J^2^^)-> the appropriate transformation O^^^ exists and, for both values of a and all nonzero 
I, has the same form 



Ki)V 



0(1) 



/ 1 











1 






1 



V2 V2 







1 

4^ 



1 







1 



1 

V2 





(75) 



On applying this transformation, >Cf take on the following form 
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V 







(76) 



where denotes the 3x3 null matrix, and f)'^'^^^ i j'^^^ are given by the formulae 

/ A - 



V 



<x(l) 





V 








S g2 

A, - (i±?^ 



(77) 



where 



Ji 
is 

J4 



(1) 



Jl 


J3 


J3 




32 




J3 


J4 


J2 



(78) 



l-2(3S-l{l + 3(3){S^ + R'') 
1 + 2/55- 1(1 + 3/5) (52 + 6i?2 

^[6/?-(l + 3/?)5] 
-i(l + 3/5)i?2 



(79) 



From (frB|) one can see that each subspace J^"^^^ splits now into two invariant subspaces 
corresponding to the following pair of triplets of new components 



^k2 



1 

JctI _ 1 



,<t/ 



,al 



^ I '^k2 '^kZ, 



(80) 



^fc3 = 75 i< + <3; 



and 



kl 
lal 

k2 

lal _ _^ 

^fc3 — 



U 



U 




'1) 



Before making another step further, let us recapitulate the operations performed so far. 
Our aim is to solve the linear equation (|Fk^ for consecutive corrections (fc = 1, 2, . . .) to 
the order parameter of the positive and negative (cr = ±1) wedge disclination line. The 
preliminary considerations have shown that the space 9 in the basis {E^, a = 1, ... ,5} 
splits into the simple sum of two subspaces Ai^^^ i Ai^'^^ invariant with respect to the 
action of therein. The Fourier analysis of that part of which belongs to Ai^^^ has 

Next, on carrying out the orthogonal 



led to the structure of the space J-'^^^ = 0^o-^i 



transformation ()75|1 in each subspace J-'[ with Z 7^ 0, we arrive at the final form of the space 
J^W which now is the simple sum comprising the subspace ^o*"^'', with Cq^^\ given by (|7I|l . 
representing the action of C'^^^^ in this subspace, and two three-dimensional nonequivalent 
subspaces for each / 7^ where the action of is represented by V^^^^ + J'^^\ see (TfTjl . 
fl78|) . The next step consists in carrying out an analogous Fourier transformation for that 
part of which belongs to the subspace A^*^^\ i.e. for components and g^g. 
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Yet, there will be no full analogy in the Fourier expansion for components {q^ay cl = 4:, 5}. 
On closer inspection of the r.h.s. terms in Eq. (jU^ . one can ascertain that the appropriate 
expansion is one comprising sines and cosines of the odd multiples of |, i.e. 



(!ka{S, 



oo 

vi:[s] + 

1=1 



Vkais) COS 



(2/ - 1)0' 



ullis) sin 



'(2/- 1)0' 



^2) 



Similarly, the module J-"*-^^ generated by that expansion splits into the simple sum of the 
subspaces invariant with respect to the action of C'^'^^^ therein. 



(2) 



1=0 



a{2) 



which have the following structure 







a(2) 



for / = 0, and 



(83) 



a(2) 



SEl cos 



SEl sin 



'(2/- 1)0' 



■(2/ - 1)0' 



SEl cos 



SEl sin 



■(2/-1) 



■(2/ - 1)0' 



^5) 



for / 7^ 0. So, the subspace JFq comprises components {f^a, a = 4, 5}, whereas in the four- 
dimensional subspace Ti'^'^'^ the components {f m^^, a = 4, 5} are to be found. The action 
of in Ti'^'^^ is represented by For / ^ the 4x4 matrices of these operators 

have the following block structure 



<x(2) 



(86) 



" (2) " (2) " (2) 

where / in the above formula denotes a 2 x 2 identity matrix, whereas J7i ,1/2 and jTa 
represent the sub matrices of J7i, J2 and jTs, respectively, comprising the last two rows and 
columns. The operator C^^'^^ has the following form 



Xa{2) 



A J - -if) + J, 



K2) 



(87) 



From (|H7j) . one can see that this operator is diagonal and again independent of o. 

Likewise, the next step to take is the partial diagonalization of with Z 7^ 0. In each 
of the subspaces Ti'^'^'' we carry out such an orthogonal transformation O'^'^'^ in order to get 
rid of the off-diagonal blocks in L"^'^^ . One can easily find that the right transformation 
O'^^^ is identical for both values of a and all / 7^ 
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0(2) 



/ 7S 7S 




V 



\ 

1 

^ ^ 
75 73 / 



(88) 



On applying this transformation to the operators Cj^'^\ they take on the form 



V 











V 



a(2) 



(89) 



where denotes the 2x2 null matrix, whereas Pj^*-^'* and are defined by the following 
formulae 



a{2) 



. _ [i-<t(2;-i)]^ 



A. 





[l+a(2/-l)]^ 
4s2 



(90) 



^ ~ 3/5/? l-f3S-l{l + 3f3){S^ + 3R'') )■ ' 

From (jHn|) one can clearly see that each subspace J^^^"^^ splits into two invariant subspaces 
corresponding to the following pairs of new components 




and 



(92) 



(93) 



This ends the first two stages of the construction, both of which were focused exclusively 
on the l.h.s. of Eqs. ()63p. The last natural stage is, therefore, to transfer all the trans- 
formations performed to the r.h.s. of those equations. The right order of carrying out all 
these transformations and the transformations themselves follow clearly from the context 
of the analysis presented above. So, the only task remaining is to introduce appropriate 
and, we hope, clear notation. Vector [A/'^J{a=i,2,3} belongs to J^^^\ whereas [A/'^J {0=4,5} - 
to AA'^'^\ The Fourier analysis in these spaces leads to the spaces !F^^^ and JF^^) where the 
components V^a^(s) i U^l{s) (a = 1, 2, . . . , 5) of the above vectors belong. They are given by 
the following formulae (see jZH), {721), (IH1|), (IH3)) 



Vk^is) = I tML{s. <P) cos{l<P)d<P, (94) 
t^.1(^) = i/o"Ar^a(^,0)sin(/0)# 



for a = 1, 2, 3, and 
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V,^1{S) = i cms, <P) COS ((^) #, 



(95) 



for a = 4, 5, and the foUowing inclusions hold 

^k{l) = [^fca°] {1=1,2,3} G -^0^^\ 
^fe(l) = ^fci]{a=l,2,3} G 



/\ro-0 
Jal 



[Vka]{a=AM ^ ^0 



a{2) 



(96) 



^fc(2) = [^fea5 f^fea]{a=4,5} ^ JF; 



a(2) 



The last step consists in taking the orthogonal transformations O^^^ and O^"^^ {I ^ 0) into 
account. It is clear that the right transformation rules for N^^■^ and N^2) have the form 



Kt) = io^'^YN^kU- 

In conclusion, the solution of the original Eq. in each order k = 1,2,.. 

perturbative expansion has been reduced to the solution of the following, in general, 
set of ordinary differential equations 



(97) 

of the 
infinite 



/^o"(1)^(tO . 

f^la{2) 
I 'ik{2) 



^^k{2)-> 
■ ^^k{2) 



dla / ^ 0, 
dla / ^ 0, 



(98) 



where 



^k(2) = bfca]{a=4,5}, 
Ika = [%a'%aJ{a=4,5} 



[Vka]{a=l,2,3}, 

Kl,<%a=i,2,3} dla 1^0, 



(99) 



dla / ^ 0, 

Fortunately, only a finite number of the inhomogeneous terms in ()98|1 are non-zero. For the 
detailed discussion of the first order corrections to the order parameter for both the positive 
and negative wedge disclination lines see Appendix 17. II 

Now the construction of the solutions for the remaining four cases ()53II56|1 ) representing 
the twist disclination hues is yet to be discussed. Let denote the solution corresponding 
to these cases. In the zeroth order of the perturbative expansion the order parameter can 
be written as follows 



where (0, \I^) is one of the following pairs 



TT 

2'« 



TT TT 

2' 4 



TT 71 

2' 2 



2'Tj 



(100) 



(101) 



So, the relevant equations for the consecutive corrections to the order parameter have 
the form 



17 



dAqi + Iqi + 3(3 {q^qi + q^q^) - \(\ + 3/3)Tr (q^qi) q^- 
-\{\ + 3/3)Tr (q^)\i - 2/3Tr (q^qi) I = J^j^. 

rp 

Let us multiply either side of the above equation by ((9i(6,^)) and (9i(9,\E') from the 
left and right side, respectively. Then, one obtains 

~dk~dkqt + Iqt + 3/? [qfqt + qUi) " W + 3/3)Tr [qfqt) qt' 



(103) 



k 1 



where 



I, = (Oi(e,vl/))^g,^Oi(e,vl/), (104) 

= me,^)fAf,^o,{e,^), (105) 

or in symbolic notation 

C'-lt=kt- (106) 
As a result, we have reduced the problem of solving Eqs. p02p for q^ to the previously 

discussed problem of solving Eqs. ()106|1 for the modified corrections % , with exactly the 
same linear operator as for the positive wedge disclination line. So, the main part of 
the analysis will proceed as it did for the case of a = +1, the only difference lying in the 
appropriately adjusted notation; all quantities corresponding to the corrections to the order 
parameter should be fitted with a tilde and superscribed with ± instead of a. For the first 
order correction to the order parameter of the twist disclination line with (B, \I') = (f ,0) 
see Appendix 17.21 

The general solution for the corrections qk {k = 1,2,...) in each of the cases discussed 
above consists of two terms, i.e. of a particular solution of the relevant inhomogeneous 
equation and a general solution of its auxiliary homogeneous partner, the latter spanning 
an infinite-dimensional space. That space splits into the subspaces corresponding each to 
the fixed value of / = 0,2,... and one of the indices (1) or (2). Obviously, only those 
elements which are continuous at s = and finite for s — oo are physically relevant and 
should be incorporated in the solution. In general, they correspond to different forms of 
excitations of the disclination line, the zero modes fj {i = 1,2,3), 9^ (a = 1,2,3) among 
them enjoying a particularly simple physical interpretation. In the first order we set the 
contribution from these solutions to zero. One can check, not easily though, that this is 
consistent with the constraints ()4H1 for k = 2. So, the final solution for the rectilinear 
disclination line up to the first order of the expansion has the form 

f = q^ + eqf, (107) 

where S G {+, — , is to refer to the different types of the disclination line, see Appendix. 
The solution p07j) may be called the minimal solution in the sense that it contains no 
excitations except for those that are forced by the anisotropic perturbation. 

For the sake of completness of the investigation, we have found numerical solutions for 
the first order corrections to the order parameter of the wedge disclination lines in the 
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special high symmetry case of /3 = 3. Then, S = 1 and the only non- vanishing functions 
P;^(s) (see Appendix 17. 1|1 satisfy the following equations 

pr + ^pr - ^P3^ - (1 + R')p,^ = -I [r" + ^R' - ^i?) , (108) 

Ps" + -/s' - ^,Ps - (1 + R')Ps = -I {r" - Ir' + • (109) 

The numerical solutions were found with the help of the Mathematica 3. package. Below, 
the zeroth order structure function R{s) as well as the first order structure functions P^^is) 
together with their first derivatives are depicted. 



1 




Figure 1: Functions R{s) and R'{s). The numerical solution has been generated within the 
interval (0.0001, 18); for s G (0, 0.0001) the numerical solution has been appended with the 
polynomial approximation smoothly matched at s = 0.0001. 



0.01 




-0.03 ' ' ' ' I 

3 6 9 12 15 

S 



Figure 2: Functions P^'is) and Pr^'{s). The numerical solution has been found in the interval 
(0.2,15); for S G (0,0.2) the numerical solution has been appended with the polynomial 
approximation smoothly matched at s = 0.2. 
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0.01 




Figure 3: Functions P^is) and P^'{s). The numerical solution has been generated in 
the interval (0.5,18); for s G (0,0.5) the numerical solution has been appended with the 
polynomial approximation smoothly matched at s = 0.5. 

The graphs above show that the first order corrections are of order ~ O.Ole. Taking into 
account the estimated values of the perturbative parameter (see the beginnig of this section) 
and the order of magnitude of the unperturbed contribution, one can see that the ratio of the 
first to zeroth order contribution is ~ 0.01. Therefore, the idea of the perturbative scheme 
finds quite a good justification. Admittedly, the above conclusion holds for the special case 
of /5 = |, but the equations for the functions P{s) giving the first order corrections to the 
order parameter (see Appendix) depend smoothly on f3 and, therefore, one can expect this 
conclusion to hold even in the generic case when /? 7^ ^• 



4 Corrections to the free energy 

In general, the knowledge of the corrections to the order parameter up to the k-th order of 
the expansion in the discussed model allows one to calculate the [k + l)-th order correction 
to the free energy. The free energy of the smooth disclination line contained in the concentric 
cylinder of the radius L^o per unit length is given by the formula 



~ 9 



(110) 



where 



Kb 



dkqijdkqij + e (piqikdjqjk — ^dkQijdkqij 



Tr (g2) + 3/3Tr (f) - ^(1 + 3/3) (Tr {q^)f] } 



1111 



is a dimensionless quantity, and the subscript 'b' is meant to refer to the smoothness of the 
underlying configuration. Inserting the expansion (j2l|l into ()111|1 . we obtain 



Kb 



KbO + SKbl +e Kb2 + ... 



;ii2) 



where 
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/«6fc = 77" / d'^sfbk (113) 

for k = 1,2, fbk represent the dimensionless corrections to the free energy density in the 
consecutive orders of the expansion. In the first two orders they have the form 



fbO 
fbl 



-dkQoijdkQoij 



Tr(g2)+3/5Tr( 



16 



m (Tr { 



Ip^ 



'^dkqoijdkqiij + diq^ikd^q^jk — -dkqoijdkqoij 



2Tr {qoqi) + 9/3Tr (f,q,) - -(1 + 3/3)Tr (g^) Tr (mi] 



(114) 



(115) 



The zeroth order contribution had been thoroughly discussed in (TUj. Let us quote here 
the most relevant results from pUj. The zeroth order contribution to Kb can be written in 
the form 



rso rL rL 

Kbo = / dssfbois) + / dssf^Q^is) + / dss 

Jo J Sq J Sq 



where 



where has the following form 



4-1(1 



fbo{s) 



9a?7o 



2*^ ' 



;ii6) 



;ii7) 



9 



aril 



'S" + 3R'' + ^R'--S' 
^ 3 



s 

— 2R^ o«c/qd2 c2^ I ^ /^ , oP\/c2 , o r>2\2 



2I3S ( - ^S'j + —(1 + 3l3)iS' + 3R' 



118) 



The regularized free energy density is given by the formula 



fbo'is) = (Sf + 3 (R'f - 3ii^M) + 1(1 + 3^)(i _ i?)2(i + ^)2+ 

+ ^(1 - Sy [32/5 + 8/3(1 + 5) + (1 + 3/3){l + Sf] - (119) 
- i(l - 5)(1 - R){1 + i?) [9/3 - 1 - (1 + 3/3)5] . 

The choice of Sq G (0, L) in pi6|) is completely arbitrary because it does not affect the value 
of Kbo- It is apparent from ()116p that the free energy of the unperturbed disclination line 
comprises two divergent terms, see the third integral in pifjj) . The first one is logarithmically 
divergent and its appearance is due to the existence of the rotational zero modes. The other 
divergent term originates from the fact that the free energy density of the ordered phase 
takes a negative value relative to that of the isotropic phase which has been set to zero. 

The asymptotic behaviour of the free energy density of the singular disclination line in 
the zeroth order has the following form 

fco{s) = ^-^{l + f3), (120) 

S o 

where the subscript 'c' had been meant to refer to singular configurations. Hence, the 
corresponding zeroth order contribution to Hc is given by the formula 
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Kco = f dssfcois), (121) 

where Sco, representing the dimensionless radius of the singular core in the zeroth order of 
the perturbative expansion, has the form 

= — — -:■ (122) 

The form of the first order correction to the free energy can be readily obtained by 
making use of the following identity 

/ (fsdiqoijdiqkij = f d^s ^Tr {q^qk) + 3/3Tr (qlqi^j - -il + 3/?)Tr [ql] Tr {q^qt 
Jc{o,L) Jc{o,L) L3 ^ ^ 4 ^ ^ 



+ / Ti (qkdiqo) da\ (123) 

which can be derived from Eq (j^3j) by first multiplying either side of it by g^, then taking 
the trace and integrating over the area of a circle C(0,L). The first order correction can 
then be written as 

2nKti= [ d^sfUs)+3f Tt (q^diqo) da\ (124) 



/C(0,L) J9C(0,L) 

where 



fbi{s) = ^ 



diqoadjqoji — -diqoijdiqoij 



(125) 



The formula ()124j) shows that to obtain Ubi one does not need to know the first order 
correction to the order parameter in detail. The only relevant information concerns the 
asymptotic behaviour of qi, for qf and q* discussed in Appendix 17.11 and 17. 2| respectively, 
the surface contribution in ()124|) vanishes at least as in the limit L oo. Now, taking 
advantage of the polynomial approximation for the functions R,S developed in [lOj, we 
can estimate the first order corrections to the free energy of the smooth disclination lines 
discussed in Appendices. Using (jl24|l . (jl25|l . and the results from [TO], one obtains 



<i = -l(w„-ir + 0^ + 0(1-'), (126) 

Firstly, from ()126I127|) it follows clearly that the perturbation removes the degeneracy. Up 
to the first order the free energy gap between the positive and negative wedge disclination 
lines is given by 

4i - Hi = 3 + 0{L-^) > 0. (128) 

But it is the twist disclination line that represents the configuration of the lowest free energy. 
This is due to the logarithmic divergence of the corresponding free energy density. That 
conclusion is in perfect agreement with the result obtained in p2] , where the corresponding 
solutions within the framework of the Oseen-Zocher-Frank theory had been considered. 
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The first order corrections to the free energy of the corresponding singular disclination 
hnes one can find remembering that in that case the only 'active' degrees of freedom are 
those connected with the average orientation of nematogenic molecules, i.e. the director 
field, whereas the degree of the orientational ordering itself is fixed at the vacuum value. 
Formally, the transition from the smooth configuration to the singular one is equivalent to 
taking the limit i?, S" — > 1. For both the wedge disclination lines this leads immediately to 
the following result 

which gives 

ACi(/?) = <i - <i = 'llMP) - I? + 0{L-') (130) 

Due to the presence of the logarithmically divergent term in ()127j) . the limit R, S —>■ 1 
cannot be carried out in such a straightforward manner. In that case, one has to perform 
the transition at the level of the corresponding free energy density. Some strenuous algebra 
gives 

rAs A) = -~il + cos2<l>). (131) 
4 

for the first correction to the free energy density of the singular twist disclination line. Like 
in JO] , we assume that the boundary of the core of the singular disclination line represents 
the set of points at which the free energy density vanishes, i. e it is equal to the free energy 
density of the isotropic phase. So, the radius of the boundary as the function of (p is given 
by 

/,*(s:(0),0) = O, (132) 

where /* (s, (p) is the full free energy density of the twist disclination line. A simple calcu- 
lation gives 

s*M = Sco - les,o{l + cos2</.) + 0{e^). (133) 

o 

The first order correction be found from the obvious formal recipe 

(134) 



as \27r JO Jsci4>;e) 



e=0 



Using ()T2njl . ()T3Tjl . and ()T33jl . one obtains 

4 \scoJ 

Hence, 



<, = -^lnf— V (135) 



<M - <i - <i = - Y4^Mf3) -ir-l In + 0{L-') (136) 

In the critical value of (3, below which the smooth disclination line is energetically 
more favourable than the singular one, has been estimated to be equal to (3ct = 0.12. In 
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the present context of the perturbative expansion that value represents the solution of the 
following equation 

A^bcoiM = 0, (137) 

where A/tbco(/5) is the zeroth order contribution to the full difference between the free energy 
of the smooth disclination line and that of the singular one AK,f^{(3), where the superscript 
refers to the three different types of the disclination line, i.e. S G {+, — The general 
equation defining has the following form 

AKfM) = 0, (138) 
where Anf^^P) as well as should be replaced by their perturbative expansions in e 



kr 



A4o(/3) + sAkIM + e'AnlM + ... 



(139) 
(140) 



where, in accord with ()137|) . (3^ = = 0.12. From the first order of the expansion of 
Eq. ()138|) one obtains the formula for the first order correction to 



( dAng^ 



\ dp 



(141) 



The value of the derivative appearing within the parentheses can be easily estimated with 
the help of the results obtained in see Fig. 4. therein. The rough estimation gives 



dp 



5.1. 



(142) 



Using that estimation and the relevant results from Subsection 3.3 in ^Oj, one obtains 



Pt ~ 0.006, 

p^ ^ 0.006, (143) 
PI ^ 0.02. 

For each type of the disclination line the first order correction to is positive, which means 
that in each case the perturbation causes the stability area of the smooth disclination line to 
expand if L2 > or to shrink otherwise. Yet, even for the twist disclination line in MBBA 
the relevant correction is too weak to stabilize the smooth configuration; the corrected value 
of P*j. still does not exceed the value Pm = 0.174 which estimates the nematic-crystalline 
state boundary. 
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5 Final remarks 



Let us recapitulate the main results of the paper. 

1. The perturbative scheme presented here represents, to the best of our knowledge, the 
first analytical attempt to go beyond the L2 = approximation for the disclination 
line in the Landau-de Gennes model. In ^ only the numerical analysis of the problem 
has been carried out. The idea of the scheme rests upon the crucial assumption that 
the relevant solution can be written as a power series in the perturbative parameter 
e. We are not able to give a proof of convergence of this series. We have obtained 
the infinite system of equations, of which only the zeroth order one is nonlinear. 
The linear operator appearing in the equations for the consecutive corrections has 
exactly the same form once the zeroth order solution has been chosen. An essential 
part in the scheme play the zero modes whose existence is due to the symmetries 
of the unperturbed system. Mathematically, they represent eigenf unctions of the 
aforementioned linear operator belonging to the eigenvalue zero. Via consistency 
conditions they have been shown to select those solutions of the unperturbed system 
which are 'matching' the perturbation. 

2. With the help of the Fourier analysis of the corrections the problem of solving the 
rather complex linear partial differential equations ()28|1 has been reduced to that of 
solving a system of linear ordinary differential equations (jHEl)- The fact that just in the 
first order the number of P-functions (see Appendix) for, say, the twist disclination 
line amounts to 15, and will grow as one passes on to the higher orders, renders the 
scheme quite impractical in the second and higher orders. Yet, it has to be stressed 
that it does by no means represent any flaw of the scheme but is a pure consequence of 
the nontrivial structure of the order parameter for a disclination line in the Landau- 
de Gennes model. Moreover, one has to remember that the first order correction to 
the order parameter field is all one needs to calculate the second order correction to 
the free energy, which in view of the fact the the model per se is not exact provides 
a satisfactory approximation. Another important aspect of the scheme is that it 
represents an approach complementary with purely numerical calculations |H]. 

3. The scheme posseses yet another advantage. It seems that it can be applied to higher 
order terms in the elastic free energy density expansion. What is important, they 
would not affect the form of the linear opeartor C. Thus, the whole structure of the 
scheme would remain intact. Certainly, the inhomogeneous terms A4 would be af- 
fected, and probably some alterations of the set of the unperturbed solutions cosistent 
with the higher order term in question would appear. 
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7 Appendix — the first order corrections to the order 
parameter 



7.1 The wedge disclination line case 

Let us first find the form of tlie correction q^. Using 
obtain 



, (jnS), (EHnZI) for A; = 1, we 



(1) 



3^6 








^^1(1) 



/ ^ {R" - f i?' + ^R) \ 





V / 



l±a_ 
2 





I (5" - -sS') 







iV(^{) = for / ^ 1, 3, 



(144) 



iVif., = 0, 



1(2) - ^{(2) = for all / ^ 0, (145) 

where R", S" and R', S' denote the second and first derivatives with respect to s, respectively. 
For the nonvanishing components of ^1(1°, ^l{T)' ?l(2) ?l(2)' accept the following 
notation 



( D+ 



P 



P.' 



V2^12 , 

2'^12 ) 
2'^13 ; 



A" 



2^12 ) 
2^13 • 



Now, using (HBHHI), 



(HUHnSl), one obtains 



(146) 



f Pi+" + iPi+' + [I - 2/35 - i(l + 3/3)(52 + R^)\ P+ + |i?[6/3 - (1 + 3I3)S]P^ 

= I {S" + IS') , 



P. 



f + 2I3S - \{l + 3/?)(52 ^ 9^2)1 



+P[6/3 - (1 + 3/3)5]Pi+ = 0, 



(147) 



26 



pr + \pr - ^Ps^ + 



2/35-i(l + 3/3)(52 + i?2 



Ps^- 



+IR[6P - (1 + 3P)S] [Pt + P^) = -\ {R" + '-R' - j,R 



Pt" + \Pt' + 
+\R[Q(3 - 



^5 



4 p+ 



; + 2/5^ - 1(1 + 3/3)iS^ + 6P2)J P++ 
+ 3/3)^]P3+-i(l + 3/5)p2p+ = o, 



+ 2/35 - i(l + 3/5) (52 + 6P2)1 P++ 



(148) 



-iP[6/? - (1 + 3/?)5]P3+ - i(l + 3/?)p2p4^ 



-i (S" - is' 



r Pf" + ipf 



2/35 - 1(1 + 3P){S^ + P2)] p- + ip[6/? - (1 + 3/3)5]P2- 



1 ^S" + \S') , 



P. 



-// I 1 p+' 1 



-'2 „2 -< 2 



+ 2/?5- i(l + 3/?)(52 + 9p2 



P2- + 



+P[6/3 - (1 + 3/3)5]Pf = 0, 



(149) 



-2/35-|(l + 3/3)(52 + p2)JPg- 
+iP[6/3 - (1 + 3/3)5] (P4- + P5-) = -| (P" - f P' + ^P 



P4-" + IP4' - ^Pa + [! + 2/?5 - i(l + 3/3)(52 + 6P2)] P4-+ 
+iP[6/3 - (1 + 3p)S]P3 - i(l + 3P)R^P^ = 0, 



P. 



-// I 1 



P.-' 



P,- + 



,,2-' 5 



+ 2/35- i(l + 3/3)(52 + p2 



+ip[6/3 - (1 + 3/?)5]P3- - 1(1 + 3/3)p2P4 



'|(5"-i5' 



(150) 



Note that Eqs. (fTITjl and (fTI^ for Pf^(s) and P^i-s) are identical. With the help of the 
formulae (jUUHT^ . (jHUj) . p46|) one arrives at the final form of the corrections 



( Pt + Pt+ 

+ {Pt + P2) COS(/)+ 

+P5+COS20 



P2 sin 0+ 
+P5+ sin 2<; 



Pg*" sin (f)+ 
+P5+ sin 20 



P^ - Pt+ 



P2 I cos ( 



-P5+COS20 



2(Pl+ + P3+COS0j / 



, (151) 
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/ Pi + P2 COS (f)+ 

COS 20+ 



+P3 cos 30 

+P4~ COS 



v 



—P2 sin0+ 
sin 20- 
—P4 sin 40 







— sin 0+ 
+P5~ sin 20- 
— P4~ sin 40 

Pf — P2" cos 0+ 



— P5" cos; 
+P3" cos 30- 
—P4 cos 40 





2 (Pf + P3" cos 30j y 



(152) 



Our considerations would be obviously far from complete if they lacked any attention to 
the relevant boundary conditions. From the physical requirement that both the corrections 
should be continuous at s = 0, one obtains the following boundary conditions for Ptis) 
{t = l,2,...,5) 



P+(0) 
P+(0) 
P3HO) 



0, 
0, 

r+ 
C4 , 

0, 



prio) 

P2i0) 

Ps'iO) 

p^io) 
n-(o) 



Cl , 

0, 
0, 
0, 
0, 



(153) 



where cf and C4 are constants to be specified in the course of the corresponding numerical 
calculations. The boundary conditions for s — > cxd are fully determined by another under- 
standable requirement that the corrections should be finite. Namely, in the limit s —>■ 00 



functions P (s) 



1, 2, . . . , 5) can be expanded in the following series 



P.^ 



±1 
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for 



1,2, 



5. 



(154) 



In the above sum all odd powers of - have been omitted. These terms vanish because the 
relevant expansions for R{s) and S{s) contain exclusively even powers of ^, see ^0]- It man- 
ifests itself through Eqs. p47H15()j) being invariant with respect to the formal transformation 
s —>■ —s. For the sake of the argument to come, let us quote here the first two expansion 
coefficients for R and S 



Rn 



Sq 



Ro 



S, 



1+15/3 



12/3(2+3/3) • 



1-3/3 



(155) 



4/3(2+3/3) • 

Inserting ()154|) together with the appropriate series for R and S into Eqs. p47H150'j) . one 
obtains the following system of equations for the expansion coefficients P^^ 



1 + 15/3)Pf + (1 - 3/5)P2^'^ = 0, 
1 - 3/3)Pi±° + (1 + 3/3)P2±° = 0, 



3±0 



1 + 15l3)Pt'' + (1 - 3l3)Pt" + (1 - mP^ 
1 - 3/3)P3±° + (1 + 3/5)P4±° + (1 + 3/5)P5±° 
1 - 3/3)P3±° + (1 + 3/3)P4±° + (1 + 3/3)P5±°. 



(156) 
(157) 
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As is the Cramer system, we get Pi = P2 = 0. In p57j) the two last equations 



are identical, so the system does not determine the values of coefficients 
uniquelly, but it is equivalent to the following system 



ctO 



4,5) 



P, 



±0 



Pt' + 



±0 



0. 



(158) 



For unique determination of P^ 



±0 



4, 5) one has to go to the next order of the asymptotic 



expansion. Equating the coefficients multiplying on either side of Eqs. ()147H15(?|l the term 
and taking ()158|1 into account, one arrives at the following systems of linear equations 



1 + 15/3)P3+2 + (1 - ?>[3)Pr + (1 - 3/5)P. 



3 -r - ^/^j-f 4 -r l-L - ^A^;-! 5 

;i - mpt' + (1 + mpr + (1 + mp. 



+2 



1 

'2' 



+2 



2I3S2 - \{l + 3/?)(^2 + 6P2) - 2(1 + 3/3)i?2P5+°, 



+2 



+ (1 + ?,(3)Pt'' 
2/352- i(l + 3/3) (^2 + 6^2) 



;i - 3/3)P3+' + (1 + W)Pt'' + (1 + 'iP)Pt = -8Pt+ 



p+O 



2(1 + 3/3)^2^^4 



+0 



(1 + mp,-' + (1 - 3/?)P4~2 + (1 - mPn 

(1 - 3/3)P3"2 + (1 + 3/?)P4-2 + (1 + 3/?)P5- 
+2 [2/3^2- 1(1 + 3/3) (52 + 6P2 

(l-3/3)P3-2 + (l + 3/3)P4-2 + (l 
+2 



3 

2' 



-32P4-°+ 



P4-°-2(l + 3/3)P2P5"°, 



3/3)P5-2 



-8Pr°- 



(159) 



(160) 



2/352 - 1(1 + 3/3) (52 + 6P2)J - 2(1 + 3/3)^2^4"°- 

Again, the l.h.s.'s of the two last equations in ()159|) and ()160|) are identical, whereas, in 
general, their r.h.s.'s differ. Hence, for the systems to be consistent one has to demand that 
their relevant r.h.s.'s should be equal. This requirement leads to the following system 



P4+° + 3P5+° 



5P4-'^ 



p.-° 



0, 
0, 



(161) 



which together with the second equation in p58j) give uniquelyP 



±0 



P^'^ = 0. To conclude. 



the consistent boundary conditions for P^ (s) 



1,2, 



P^(oo) 



0, 



, 5) have the following form 
1,2. ..,5. (162) 



7.2 The twist disclination line case 

~ 

For the non- vanishing components of Afi , using (jS21), (IHSI), dMUHZI), p()4|l . and p()5|l with 

k = 1, one obtains 







V 







(163) 
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( 












TJlLl 




'^12 




\ *^13 


} 



( ^ 



->f±2 























V uif ) 





cos2^ (i?" + ii?' - ^i?) \ 

-^cos2* {S"-lS' 
-^sin2^^ {R" + ^R'- j,R 
-|sin2* [S" + ^^S' 
-^sin2^ (5"- i^' 





\ 




^12 






yil.2 






7-~r/±2 
'^ll 






^\2 








I 






\ 









^C0S2V1/ i^R"-lR' + ^R 

-3^ COS 2* {S"-lS') 


^R" - f i?' + ^R 
^sin2V(5"-i5'f 



(164) 



(165) 



(166) 



^/±2 
^1(2) 



^1(2) 



yiL2 
•^15 

1 rr/±2 
\ '^15 



( -\cos^ (r" + \R' - j,R 



-isin^ (R" 



S" -^S' 

s 

-^cos^(5"- 
cos^ (i?" - fi?' 
^sin^(5''-i5') 



-^sm 



(167) 



(168) 



where takes on one of the four values {0, |,|,^} (for each configuration considered 
here © = |). One should remember that in the above formulae the symbol ' denotes the 
differentiation with respect to s exclusivly for the functions R or S. Otherwise, it refers to 
the orthogonal transformations (9*^^^ and O*^^) in the submodules J-'^^^'' and ^"^^ [1^0), 
respectively. The system of equations pertinent to the discussed case has the form 



where 



Q'fc(i) 



''-o yfc(2) 

/'+(2)^^' 
%{2) 



N 



.±0 

N 



fc(2) 



= N 



for / 7^ 0, 



for I ^ 0, 



(169) 
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:?-L0 _ 
Ikil) = 

Qk{i) = 

it±0 _ 
Qk{2) = 

Ika = 



„-,-LOl 



Jka\{a=l,2,3}, 
.^fea°]{a=4,5}, 



for /^O, 



for / ^ 0, 



(170) 



for k = 1,2, . . .. Specifically, for k = 1, using ()163H17(]j) . one can obtain the final forms of 
corrections Qi and write the analogous of Eqs. ()147H15()1l for all the four cases. Due to the 
length of the resulting formulae we confine ourselves to quoting the result for one case only, 
i.e. for \I' = 0. Let us replace the superscript by * in all quantities referring to that case. 

Now, for the non-vanishmg components of Qi(2) i 1i{2) T 0) mtroduce the 

following notation 



Pi 

JD* 1 

^2 — ^ 



P* 



J_„-,*0 

j_7*i 
1 

1 

2'^13 ' 

1 y*2 



1-1*2 
2^12 5 
l~/*2 



P* — J_,~,'*3 
^9 - 76^11 ' 

p* 1^,1*3 

-'10 ~ 2 12 ' 
P* — ii~/*3 
11 ~ 2'^13 ' 
1 ~/*2 
2^14 5 
1 ~/*2 
2'^15 ' 

2^14 ' 

2'^15 • 



P* 
-'12 

P* 
^13 

P* 
-'14 

^15 



(171) 



Using ()163H17lj) . one arrives at the following systems of ordinary differential equations for 
P:is) (^ = 1,2,.. .,15) 



Pi 



-kii I 1 r>*i 



- i(l + 3/?)(S2 + P?)\ P{ + ii?[6/? - (1 + 3/3)^]P2* 
= h [S" + -sS') , 



Pf + \P^' - jjPI + [I + 2I3S - 1(1 + 3/?)(52 + 9P2)] P2*+ 
+i?[6/3 - (1 + ^I3)S]P^ = 0, 

P*" _L 1 P*' J_ P* _L 



-2/35-i(l + 3/5)(52 + i?2^ 

y/ I 1 D/ 



+ip[6/? -il + 3f3)S] {PI + P5*) = ^ ( P" + iP' 



P4*" + 7P4' + [1 + 2/35 - 1(1 + 3/3)(52 + 6P2)] P4^+ 
-iP[6/5 - (1 + 3(3)S]P^ - \{l + 3/5)p2P5* = 3^ + i^' 



JD*II I 1 p*/ 4 p* 

-^5 s-f^5 ^2-^5 



+ - i(l + 3/?)(S2 + 6P2)] P5*+ 
[ +ip[6/3 - (1 + 3/3)5]P3^ - 1(1 + mR'Pl = -k - -s^' 



(172) 



(173) 
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s 6 G I 



1-2/35- + 3/3) (52 + i?2 



+\R[Q(3 - (1 + 3(3)S] (P; + Pg' 



-1-2 iS" 



P*" + iP;' - ^P; + [I + 2/?S - 1(1 + 3/?)(52 + 6P2)] p*+ 
+ip[6/3 - (1 + 3/3)5]P6^ - i(l + 3/3)p2Pg- = 1 (^R" + ii?' _ , 



^Pg* + [f + 2/35 - 1(1 + 3/3) (^2 + 6P2)] p*+ 
+ip[6/3 - (1 + 3/3)5]P6^ - 1(1 + 3/3)p2p; = 1 (^p" _ |p' + ^p' 



-2/35-i(l + 3/3)(52 + p2 



Pq* + 



+iP[6/3 - (1 + 3/?)5] (P,*o + P,\) = (P" - f P' + 



+ -sPio - ^Pio + [i + 2/5^ - 1(1 + 3/9)(S2 + 6P2)] P,^o+ 
+ip[6/3 - (1 + 3/3)5]P9^ - i(l + 3/3)p2p,\ = (5" - 



-^11 + s-^ll 



16 p* I 

s2 Al + 



+ 2/3^ - i(l + 3/3)(^2 ^ gp2)l p. ^ 



+ip[6/3 - (1 + 3/3)5]P9^ - 1(1 + 3/3)p2Pi*o = 0, 



p*// 

-'12 



+ iPil - j,P^, + 



P*" 
-'13 ^ 



-'14 



1 p*/ 

s-^13 



+3/3PPi^3 



^2-^13 ^ 



- /?S - i(l + 3/?)(^2 ^ 3p2)] p.^^ 



- /?^ - 1(1 + 3/?) (^2 ^ 3p2)J p.^^ 

+3/3PPi*2 = ^ {S" - IS') , 



iPi1-M'4+ |-/35-|(l + 3/3)(52 + 3p2) Pi^,+ 



+3/3PPi^5 = {S" - i^' 



p*// lip 

-^15 +7-^] 



s2 -f^l5 + 



(3S - i(l + 3(3){S^ + 3P2)1 Pi*5+ 



+3/?PP^i = I (P" - f P' + ^P 



(174) 



(175) 



(176) 



(177) 



Using dill), (jnni), (EH), (|HUHH^. (|n2Hni), (Uni, and dnH), one obtains the final form of g* 
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-2 [P* + P* cos 
+P^ cos 20 + Pg* cos 



— P*2 sin 0+ 
+ (Pi^3-Pi*4)sin20+ 
+P]*5 sin 30 



P*2 sin 0+ 



+ (Pi^3-Pi^4)sin20+ 
+P*5 sin 30 



Pi" - p:+ 

+ {P^-P^_-P^) COS0+ 

+ (P(r-P5"'-Pi"o)cos20+ 

+ (P9*-P8*)COS30- 



-P]*^ COS' 



-p 



13 



12 COS0- 

P*.) COS 20- 



-p 



15 



cos; 



(P; - P*)sin0+ 
+ (Pi^o-n')sm20- 
— Pg sin 30 — P]*]^ sin 40 



— P*2 cos 0— 
■(Pi^3 + Pi*Jcos20- 

— P*5 COS 



(P7* - P2*)sin0+ 

+ (Pi^o-n*) sin 20- 
-Pg sin 30 — P]*i sin 

Pi* + P4*+ 
+ (P3* + P2* + Py*) COS0+ 
f (P6* + P5* + Pi*o)cos204 



+ 



V 



;P(r"'+P5* + Pi*o)cos20+ 
+ (P* + P*) COS 30+ 
+Pi\ COS 40 



In a way analogous to that described in detail in Appendix 17. ![ one can obtain the consistent 
boundary conditions for the functions P*(s) = 1, 2, . . . , 15). They read 



(178) 
(179) 

(180) 

(181) 

(182) 
(183) 



r pm = 1 

1 A^(o = 0, 1 


' Pi*(oo) 
. P^(oo) 


= 0, 
= 0, 


f ^3^(0) = 0, 
pm = ci 
[ pm = 0, 


r P3*(00) 

n(oo) 
I n^oo) 


= 0, 
= 0, 
= 0, 


f ^6 (0) = 0, 

< p;(o) = o, < 

[ pm = 0, 


r P6*(00) 

P;(oo) 

. i^8(00) 


= 0, 
1 

_ 16'^ 

16 


f ^9 (0) = 0, 

< Pi*o(0) = 0, < 
[ Pu{0) = 0, 


f n^(oo) 


= 0, 
= 0, 
= 0, 


/ pm = 0, 1 
1 aijW = 0, 1 


f ^1*2(00) 


1 

32 
32 



PiM 



pm 



0, 
0, 



Puioo) = 
PUoo) 



T 

32' 



where Ci and C4 are constants to be determined in the course of numerically solving Eqs. 
fl1V'il1V3|) . 
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